Majorana spin liquids and projective realization of SU(2) spin symmetry 
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We revisit the fermionic parton approach to 5 = 1/2 quantum spin liquids with SU(2) spin 
rotation symmetry, and the associated projective symmetry group (PSG) classification. We point 
out that the existing PSG classification is incomplete; upon completing it, we find spin liquid states 
with S = 1 and 5* = Majorana fermion excitations coupled to a deconfined Z2 gauge field. The 
crucial observation leading us to this result is that, like space group and time reversal symmetries, 
spin rotations can act projectively on the fermionic partons; that is, a spin rotation may be realized 
by simultaneous SU(2) spin and gauge rotations. We show that there are only two realizations of 
spin rotations acting on fermionic partons: the familiar naive realization where spin rotation is not 
accompanied by any gauge transformation, and a single type of projective realization. We discuss 
the PSG classification for states with projective spin rotations. To illustrate these results, we show 
that there are four such PSGs on the two-dimensional square lattice. We study the properties of 
the corresponding states, finding that one - with gapless Fermi points - is a stable phase beyond 
mean-field theory. In this phase, depending on parameters, a small Zeeman magnetic field can 
open a partial gap for the Majorana fermion excitations. Moreover, there are nearby gapped phases 
supporting Z2 vortex excitations obeying non-Abelian statistics. We conclude with a discussion of 
various open issues, including the challenging question of where such S = 1 Majorana spin liquids 
may occur in models and in real systems. 



I. INTRODUCTION 

Some of the most intriguing states of matter are those 
beyond the conventional classification of phases accord- 
ing to spontaneously broken symmetry, band theory, 
and Fermi liquid theory. The classic examples of such 
so-called exotic states are the fractional quantum Hall 
liquids, which are characterized by topological order^ 
and associated properties such as fractionally charged 
excitations.^ Quantum spin liquids^ - ground states of 
Mott insulators with no spontaneously broken symmetry 
- are another class of states that are in many cases ex- 
otic, and are often also characterized by topological order 
and fractionalized excitations. ^^^"^ A number of experi- 
ments over the last several years have uncovered materi- 
als where exotic quantum spin liquids may be present. ^'^^ 

In systems with S = 1/2 local moments and SU(2) spin 
rotation symmetry, S = 1/2 spinous generally play an 
important role in the theory of spin liquid states, whether 
as quasiparticle excitations, as formal objects in terms of 
which the theory is constructed, or both. Such spinous 
can be bosons or fermions, and can also obey fractional 
or non-Abelian statistics. 

Several recent works have raised the fascinating 
prospect of spin liquids where spinon quasiparticles do 
not carry S = 1/2, but are instead S = 1 Ma- 
jorana fermions. ^^"^^ In particular, such states were 
shown to occur in exactly solvable SU(2)-invariant 
generalizations^"'^'^^^^'^'"'^^ of Kitaev's honeycomb lattice 
model. Of course, such models are rather special, and 
an approach to study 5 = 1 Majorana spin liquids in 
more general 5 = 1/2 spin models is a desirable comple- 
ment to the exact solutions. Biswas et. al. have made 
a fascinating proposal in this direction, constructing a 
mean-field theory of such a state on the triangular lattice. 



We note that Majorana fermions also play an important 
role in other spin liquids of recent interest with SU(2) 
spin symmetry.^^'^^ 

In this paper we show that, surprisingly, 5 = 1 Ma- 
jorana spin liquids fit naturally into a well-known con- 
struction of spin liquids based on 5 = 1/2 fermionic par- 
tons. Such spin liquid states can be classified in terms 
of their projective symmetry group (PSG).^^ We show 
that the existing PSG classification of 5 = 1/2 fermionic 
parton states is incomplete; upon completing it, we find 
the 5 = 1 Majorana spin liquids. The spin liquids we 
find have, in addition to a triplet of 5 = 1 Majorana 
fermions, a spin singlet Majorana fermion. Going be- 
yond mean- field theory, all of these fermions are coupled 
to a deconfined Z2 gauge field. Moreover, wavefunctions 
for these states are easily obtained via Gutzwiller projec- 
tion. 

Before delving into details, let us first make the above 
assertions plausible. In the 5 = 1/2 fermionic parton 
approach, space group and time reversal symmetries act 
projectively on fermions - that is, such operations are 
realized as a product of the naive operation combined 
with an appropriate SU(2) gauge transformation.^^ The 
crucial observation, which to our knowledge has not been 
made before, is that spin rotations can also act projec- 
tively on the 5 = 1/2 partons. We show that there are 
only two distinct realizations of spin rotation symmetry 
- a familiar naive realization, where spin rotations alone 
are a symmetry, and a projective realization, where only 
combined spin and gauge rotations are a symmetry. The 
projective realization can be thought of as a "locking to- 
gether" of SU(2) spin symmetry and SU(2) gauge trans- 
formations, and is analogous to color-flavor locking in 
high-density quantum chromodynamics.^^ 

We find that 5 = 1 Majorana spin liquids occur when 
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spin rotations arc realized projectively. The spin symme- 
try of such states is not readily apparent if one works in 
terms of 5 = 1/2 partons. Manifest spin rotation invari- 
ance is recovered upon writing the S = 1/2 partons in 
terms of Majorana fermions that turn out to transform 
as singlets and triplets under projective spin rotations. 

When spin rotations arc realized projectively, the PSG 
classification needs to be re-done. We show that Majo- 
rana spin liquid PSGs - henceforth referred to as Majo- 
rana PSGs - are in one-to-one correspondence with the 
SU(2) PSGs in the existing classification {i.e. for the 
naive realization of spin rotation symmetry). In general, 
for a given lattice there are only a few SU(2) PSGs, and 
therefore there are not many Majorana PSGs, at least 
if we insist on keeping time-reversal invariancc and all 
lattice symmetries, as we do - for simplicity - in this 
paper. For example, on the square lattice there are four 
Majorana PSGs,^^ and on an anisotropic triangular lat- 
tice there are two.^^ On the perfect (isotropic) triangu- 
lar lattice we show that there are actually no Majorana 
PSGs. This occurs because the mean-field Hamiltonian 
for a Majorana spin liquid breaks time reversal (and usu- 
ally also reflection symmetry) if it has a non-bipartite 
structure of hoppings; therefore, unless one allows break- 
ing of some symmetries, many frustrated lattices are not 
expected to admit any Majorana spin liquids. On the 
other hand, if we broaden our scope to allow for breaking 
of time reversal and some lattice symmetries, Majorana 
spin liquids are certainly possible on frustrated lattices. 
In Sec. VIII, wc speculate on the implications of these 
observations for finding Majorana spin liquids in realis- 
tic models and in experiments. 

Beyond development of the general results mentioned 
above, we discuss Majorana spin liquids on the square 
lattice. The primary purpose of this discussion is to give 
a concrete illustration of our more general results, and to 
discuss the properties of some spin liquid states arising 
from our construction. There are four Majorana PSGs 
on the square lattice; for each of these we discuss the 
state with the simplest mean-field Hamiltonian (i.e. with 
only the shortest-ranged hopping allowed by symmetry). 
Three of these states have nested Majorana Fermi sur- 
faces, and wc expect these do not describe stable spin 
liquids beyond mean-field theory. One state, however, is 
characterized by gapless Fermi points, and we show that 
it is a stable phase. This state is dubbed the MBl-Dirac 
state. The behavior of this state in a Zeeman magnetic 
field is interesting; depending on parameters, a small Zee- 
man field opens a gap for some of the Majorana fermions. 
Similar effects of Zeeman field were noted in the exactly 
solvable model of Ref. 14. The Fermi points of the MBl- 
Dirac state acquire a full gap upon introducing either a 
weak columnar dimerization or a weak breaking of parity 
and time-reversal. In the presence of both these orders, 
it is possible to have a state where Z2 vortices are bound 
to an odd number of Majorana zero modes and carry 
non-Abelian statistics, as in the B-phase of the Kitaev 
honeycomb lattice model. 



Our mean-field theory is distinct from that of Biswas 
et. al. (Rcf. 13): in their formalism there is only a triplet 
of Majorana fermions, where in ours there is also a spin 
singlet Majorana fermion. Understanding the relation- 
ship, if any, between these two mean-field theories is an 
open problem. We contrast the two approaches in Ap- 
pendix D. Briefly, our formalism can incorporate fluctu- 
ations about mean-field theory using standard ideas of 
slave particle gauge theories, while Ref. 13 seems to re- 
quire a more novel approach, which would be interesting 
to study in detail. Projected wavefunctions can be easily 
obtained in our formalism, while it is not yet clear how 
to do this following Biswas et. al. These distinctions 
notwithstanding, it should be emphasized that the states 
obtained via these two approaches are similar in their 
physical properties, and in some cases the two mean- 
field theories may even describe two diff'erent limits of 
the same phase. 



Wc now outline the remainder of the paper. In Sec. II, 
we briefiy review the fermionic parton approach to S = 
1/2 spin liquids, and the classification of such states by 
projective symmetry group. A variety of useful notation 
is also introduced there. In Sec. Ill, we discuss the pos- 
sibility of projective spin rotation symmetry, and show 
that there are only two possible realizations of spin rota- 
tion symmetry in the fermionic parton approach, subject 
only to some minimal assumptions. Next, in Sec. IV, 
wc study the most general mean-field fermion Hamilto- 
nian invariant under projective spin rotations, and show 
that its single-particle excitations are S = I and S = 
Majorana fermions. We also discuss the low-energy ef- 
fective Z2 gauge theory of Majorana spin liquids, and 
make some comments about their projected wavefunc- 
tions. Section V is concerned with the classification of 
Majorana PSGs. In Sec. VA we establish the one-to- 
one correspondence between Majorana PSGs and SU(2) 
PSGs in the existing PSG classification. In Sec. VB we 
enumerate the four Majorana PSGs on the square lattice, 
and in Sec. V C we discuss frustrated mean-field ansatze 
for Majorana spin liquids and time reversal symmetry 
breaking. In Sec. VI, we study Majorana spin liquids 
on the square lattice at the mean-field level, considering 
each of the four PSGs. The properties of the stable MBl- 
Dirac state are considered in more detail in Sec. VII. In 
particular, we show that the MBl-Dirac state is a stable 
phase and discuss the effects of Zeeman magnetic field. 
We also consider the properties of some nearby gapped 
phases, including some with non-Abelian statistics of Z2 
vortices. Various technical details, as well as a discussion 
of the approach of Ref. 13, are contained in the appen- 
dices. The paper concludes with a discussion in Sec. VIII. 
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II. REVIEW OF FERMIONIC PARTON 
APPROACH TO S = 1/2 SPIN LIQUIDS 

We consider a system of S' = 1/2 spins placed on the 
sites r of some regular lattice, with Hamiltonian 



H 



E 

(r,r') 



(1) 



Here the sum is over distinct pairs of sites (r,r'); for 
later applications we take each such pair to be ordered 
according to some arbitrary convention. The operator 
Sj. generates rotations of the spin at site r, and satis- 
fies the commutation relations [5*^,5^,] = iS^r'^^^^ S^. 
Our focus is on systems obeying full SU(2) spin rotation 
symmetry, as well as time reversal and space group sym- 
metry. The precise form of the Hamiltonian will be less 
important for us than its symmetry, since wc arc primar- 
ily concerned with constructing and classifying possible 
states - the much more difficult problem of finding spe- 
cific models that realize the new states we identify will 
be left for future work. Even so, it is worth noting that 
in most cases of interest the exchange couplings Jrr' will 
be predominantly antiferromagnetic (positive), although 
it is not necessary for all the exchange couplings to be 
positive. Moreover, additional multi-spin exchange terms 
(not written) can also be included in the Hamiltonian. 

We are interested in constructing possible spin liquid 
ground states, which are simply ground states that pre- 
serve all of the microscopic symmetries of the original 
model. (Occasionally we will also use a looser defini- 
tion of spin liquid that only requires the preservation 
of translation and spin rotation symmetry.) One ma- 
jor approach to constructing spin liquid states begins 
by rewriting the spin operator Sr as a bilinear of par- 
tons (bosons or fermions) and ends in the construction 
of wavefunctions, as well as associated low-energy effec- 
tive field theories. For reasons discussed at the end of 
Sec. Ill, in this paper wc shall confine our attention to 
fermionic partons. In the remainder of this section we 
review the fermionic parton approach, mostly following 
Ref. 19, before proceeding to our results in Sec. III. Our 
intent is not to provide a complete review, but rather to 
remind the reader of the basic facts, emphasizing those 
aspects important for connecting to the remainder of the 
paper. To this end, we use notation difi'ering from most 
treatments in the literature. 

The spin operator is written as 



r/3' 



(2) 



where fl^ creates a spin- 1/2 fermion of spin a =t,-|. at 
site r. The fermions obey canonical anticommutation re- 
lations, <T = {a^,a'^,a^) is a vector of the 2x2 Pauli 
matrices, and summation over repeated indices is im- 
plied. In order to have a faithful representation of the 
spin model, we must also impose the local constraint 



It is clear from Eqs. (2) and (3) that there is a lo- 
cal redundancy under U(l) gauge transformations /^q, — > 
e*'^'" fra ■ In fact the full local redundancy is known to be 
SU(2).^^'^^ To expose this, we introduce the 2x2 matrix 



Fr 




(4) 



The reason for introducing this matrix, rather than work- 
ing with a two-component spinor as is more common, is 
that left-SU(2) rotations of Fr are spin rotations, while 
right-SU(2) rotations are gauge transformations. There- 
fore, when we work with objects built from ^V, transfor- 
mation properties under both spin and gauge rotations 
are manifest. The spin operator can be written 



(5) 



which is manifestly invariant under the SU(2) gauge 
transformation FrUr, where Ur € SU(2). More- 

over, gauge transformations are generated by 



1 



triFrCTF^). 



(6) 



The constraint Eq. (3) can be expressed as = 0, and 
this automatically implies also 0^ = 0^= 0. Therefore, 
we have the local SU(2) constraint 



Gr = 0. 



(7) 



A fact that will be useful later on is that Sr and Gr form 

a complete set of fermion bilincars for the single site r. 
We have the commutation relations 



[Sr, Gl. 



= 0. 



(8) 
(9) 
(10) 



To proceed, one writes down a quadratic mean-field 
Hamiltonian in terms of the fermions. This Hamilto- 
nian can be obtained by decoupling the quartic spin ex- 
change interaction, but this is not necessary and we will 
not frame our discussion in that language. The most gen- 
eral quadratic Hamiltonian satisfying [Hq, S*'] = 0, where 

5* = Sr, is 



Ho 



,ir{FrFl,) + ^rr' ^<Fr^' 4') 



(11) 



where a^ir), Xrr' ^^'^ xl-r' = 1^2,3) are real param- 
eters. A choice of these parameters is referred to as a 
mean- field ansatz. 

A spin liquid wavefunction is generated from this 
Hamiltonian by Gutzwiller projection: 



f 

J raJ r 



1. 



(3) 



(12) 
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where |V'o) is the ground state of Hq, and V implements 
GutzwiUer projection (i.e. projection onto the subspace 
satisfying the constraint Gr = 0). More precisely, V = 
llr^^' where 

From this form it is clear that Gl.\tp) = 0. 

The form of Hq guarantees that \tp) is a spin singlet: 
S'ltjjo) = 0, and therefore S^\tP) = S^V\tpa) = VSyja) = 
0.-^^ For \tp) to describe a spin liquid, it must also preserve 
time reversal and space group symmetries (or a subgroup 
of these symmetries, if wc use a looser definition of spin 
liquid). This occurs exactly when Hq is invariant un- 
der projective symmetry group (PSG) transformations.^^ 
Consider a space group operation S : r — ^ S{r). This op- 
eration acts on spin operators by 

S:Sr^Ss(r), (14) 

and we require this operation to leave \Tp) invariant (pos- 
sibly up to multiplication by an overall phase). Acting on 
fermion operators, we have the projective transformation 

S:Fr^ Fs^r)U^, (15) 

where G SU(2) is an arbitrary gauge transformation, 

which docs not affect the transformation of the gauge- 
invariant spin operators. In order for S" to be a symmetry, 
we require that there exists some choice of such that 
Eq. (15) leaves Hq invariant. Similarly, time reversal is 
implemented as an anti-unitary operation sending 

r-.Fr^ {ia^)FrU^, (16) 

where again Uj^ must be chosen to leave Hq invariant. 
Spin rotations are realized by Fr — > UFr for U G SU(2), 
and do not require any gauge transformation (but see 
Sec. III). Moreover, Hq is also invariant under a sub- 
group of pure gauge transformations, called the invariant 
gauge group (IGG). The IGG can be SU(2), U(l), Z2, or 
products of these groups. There is thus some freedom in 
the choice of and Uj^, since these gauge transforma- 
tions can always be multiplied by an element of the IGG. 
This situation is expressed by writing 

SG = PSG/IGG, (17) 

where SG stands for the symmetry group of the spin 
model. 

To summarize this discussion, and give a more pre- 
cise definition of PSG, we say that a PSG is specified by 
the following set of transformations (and also products 
of these transformations): 

S : Fr ^ Fs^r)U^ (18) 

T : Fr ^ iia^)Fr.Uj (19) 

Spin rotation : Fr UFr (20) 

IGG : Fr FrU^. (21) 



Here, S runs over all space group operations, and a 
parametrizes the elements of the IGG. We also require 
that there exists a (non-zero) ansatz invariant under 
these transformations. Moreover, we require that the 
ansatz is invariant only under those pure gauge trans- 
formations in the IGG. Two collections [Eqs. (18 - 21)] 
of such transformations are equivalent, and are consid- 
ered realizations of the same PSG, if they arc gauge- 
equivalent. (Wc do not require the ansatze associ- 
ated with two sets of transformation laws to be gauge- 
equivalent.) 

PSGs can be classified following Rcf. 19, where the 
classification was worked out for the square lattice for 
the cases IGG = SU(2),Z2, where there are only a fi- 
nite number of PSGs. The classification was also par- 
tially worked out for IGG — U(l), where there arc an 
infinite number of PSGs. PSGs are often referred to by 
type of IGG; for example, if IGG = Z2, we say that the 
PSG is a Z2 PSG. It should be remarked that a very 
similar PSG classification exists in the bosonic parton 
approach.^^ The main difference is that the gauge struc- 
ture for bosonic partons is only U(l) and not SU(2). 

So far we described how to generate a wavefunction 
from a mean-field ansatz, but we have said nothing about 
how to arrive at a low-energy effective theory. We sketch 
a prescription here, which builds on ideas introduced in 
Ref. 26, and has subsequently been used in many works. 
Beyond the brief discussion here, we also illustrate this 
prescription in greater detail, via a concrete example in 
Sec. IV. One introduces a dynamical lattice gauge field 
with gauge group given by IGG, and couples it to the 
fermions.^'' The Gauss' law constraint is chosen so that 
the gauge theory reduces to a 5 = 1/2 Heisenberg model 
in the strong coupling limit. Barring accidental fine- 
tuning - which is anyway easily corrected the resulting 
low-energy theory has precisely the global symmetries of 
the microscopic spin model of interest. Moreover it re- 
duces to a spin model in the same universality class {i.e. 
with 5 = 1/2 spins, short-range interactions, and the 
same symmetries), in the strong coupling limit. There- 
fore it is expected to be a legitimate low-energy effective 
theory, in the sense that its phases and phase transitions 
occur for some spin model in the same universality class 
as the microscopic spin model of interest. Furthermore, 
such gauge theories also arise naturally upon studying 
fluctuations about mean-field theory. 

The spin liquid phase is the deconfined phase of the 
low-energy effective gauge theory. It should be noted 
that when IGG — U(1),SU(2), there may not be a sta- 
ble deconfined phase; in that case the low-energy the- 
ory does not describe a stable spin liquid phase. When 
IGG = Z2, deconfinement of the gauge field is robust, as 
it is protected by the non-zero energy gap to Z2 vortex 
excitations. 

It is important to keep in mind that the classification 
of PSGs is not the same as a classification of spin liquid 

phases. For instance, there can be distinct spin liquids 
with the same PSG. This occurs, for example, when some 
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parameter of a mean-field ansatz can be tuned to trans- 
form a state with a fermion gap into a gapless state. 

Finally, we briefly comment on the relationship be- 
tween the projected wavefunction and effective theory 
obtained from the same mean-field ansatz. In our opin- 
ion, this issue is poorly understood and in need of more 
attention in future work. The effective theory, by de- 
sign, correctly captures the universal long-wavelength 
physics of a given spin liquid phase; there is no guar- 
antee that the wavefunction does the same. Indeed, in 
a number of cases there is compelling evidence that the 
long- wavelength properties of projected wavefunctions do 
not match the corresponding effective theory. ^^"^^ On the 
other hand, a class of wavefunctions for Z2 spin liquids 
does have the same Zi topological order as the effective 
gauge theory.^^'^^ Moreover, the wavefunctions provide 
short-distance information - energetic information, for 
example - that is inaccessible using the effective the- 
ory approach. Therefore it would be desirable to bet- 
ter understand the circumstances under which projected 
wavefunctions capture the correct long-distance behav- 
ior of the low-energy effective theory, and, in other cir- 
cumstances, to learn how the wavefunctions may be im- 
proved. 

III. PROJECTIVE REALIZATION OF SU(2) 
SPIN SYMMETRY 

Here, within the framework of the S = 1/2 fermionic 
parton approach reviewed above, we consider the possi- 
bility that spin rotations are realized projectively. That 
is, we consider the possibility that the mean-field Hamil- 
tonian is not invariant under the naive spin rotation 
Eq. (20), but is invariant when naive spin rotation is 
combined with an appropriate gauge transformation. Be- 
cause spin rotations are a continuous symmetry, the con- 
ditions on how it may be realized arc quite restrictive. We 
will show that, subject to minimal assumptions, there are 
only two distinct ways to realize spin rotation symmetry: 
the naive transformation of Eq. (20), and the projective 
spin rotation 

Spin rotation : Fr IJFrU''. (22) 

When IGG = SU(2), Eq. (20) and Eq. (22) are not 
distinct. However, in Sec. IV below, we will consider 
the most general mean-field Hamiltonian invariant un- 
der projective spin rotations and see that generically 
IGG = Z2. 

We assume spin rotations are generated by the Hermi- 
tian operators T*, and that [Ho,T''] = 0, where Hq is the 
mean-field Hamiltonian. Note that we do not assume Hq 
has the form given in Eq. (11). It follows from Noethcr's 
theorem and the fact that the mean-field Hamiltonian is 
quadratic that T' is a fermion bilinear. We make the 
following further assumptions: 

1. = ^^t;. 



2. For a gauge- invariant state (which satisfies 

G;|V) = o),r'|V) = s*|V). 

While it might conceivably be possible to relax some of 
these assumptions, we have not found a sensible way to 
do this, and we will not consider this possibility here. 

Here and throughout this paper, we will restrict atten- 
tion to mean-field ansatze that fully connect the lattice. 
This means that any two sites ri and r2 are joined by 
a path of lattice bonds (r,r') such that, for each bond 
in the path, a fermion bilinear coupling Fr with Fr' ap- 
pears in Hq with nonzero coefficient. One reason for this 
restriction is just simplicity. A deeper reason is that, 
for a mean-field ansatz where the lattice breaks into two 
or more disconnected components, the PSG classifica- 
tion reduces to a separate PSG classification for each 
disconnected component, and the IGG will be a prod- 
uct of IGGs for each of the disconnected components. 
Therefore the more basic problem is to classify PSGs 
(and construct corresponding spin liquid states) for fully 
connected mean-field ansatz. 

Since S^. and G\. are a complete set of single-site 
fermion bilinears, the most general form of satisfy- 
ing assumption (2) is 

Ti = Si. + M'jGi, (23) 

where M*-' is an arbitrary real 3x3 matrix. For assump- 
tion (3) to hold we must also have [Tl^Ti.] = ie'^'^T^; it 
is shown in Appendix A that this implies either Mj. = 
or Mr S SO (3). We can therefore make a gauge transfor- 
mation so that on every site either Mr = or M^^ = 6^^ . 
Suppose that on one site r, Mr = 0, while on another r', 
M^^i = 5^^ . In this case, there is no spin rotation invariant 
fermion bilinear coupling Fr and Fr' . Since we assume 
the ansatz is fully connected, this means we must either 
have Mr = everywhere, or M^P = 6^^ everywhere. 

We have therefore shown that the only two possibilities 
for the generator of spin rotations are T' = 5', or 

T' = S' + G\ (24) 

where 

ff=^G;. (25) 

r 

This form generates the projective spin rotations of 
Eq. (22). We shall now proceed to study mean-field 
Hamiltonians, and the corresponding spin liquid states, 
where spin rotations are realized projectively in this fash- 
ion. 

At this point, it is natural to ask whether any anal- 
ogous results hold for bosonic partons. Obviously the 
SU(2) gauge structure of fermionic partons is the crucial 
element in the above discussion, because it allows for a 
natural association, expressed in Eq. (22), of a gauge ro- 
tation with a given spin rotation. Bosonic partons have 
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only a U(l) gauge structure, so we expect that spin ro- 
tations cannot be realized projectively with bosonic par- 
tons, as long as the spin symmetry is SU(2). On the other 
hand, if the spin symmetry is only U(l), we expect that 
projective spin symmetry can be realized with bosonic 
partons. Study of this possibility is left for future work. 



IV. MAJORANA SPIN LIQUIDS 

Here, we consider the most general mean-field Hamilto- 
nian invariant under projective spin rotation symmetry 
generated by T* = S'' + G' , and show that its single- 
particle excitations are 5 = 1 and 5 = Majorana 
fermions. Moreover, we show that IGG = Z2, and 
write the low-energy effective Z2 gauge theory describing 
the spin liquid state. We also make some comments on 
projected wavefunctions. The classification of Majorana 
PSGs, and specific examples of Majorana spin liquids, 
are discussed in later sections. 

The most general quadratic Hamiltonian invariant un- 
der projective spin rotation symmetry is 

^^0= ^ [ixlr'HFr.Fl')+^xlr'H<y'F,<y^Fl')\ (26) 

(r,r') 

The sum is over ordered pairs of lattice sites (r, r'); the 
ordering of the pairs is fixed but arbitrary. If x^^, — 
for all bonds, then we have IGG = SU(2), in which case 
there is no distinction between projective and naive spin 
rotations. Therefore we always want to consider x^^, 7^ 
for some bonds. The ground state of Hq is a singlet 
under projective spin rotations, that is T'lV'o) = 0, and 
also (V'olT^IV'o) = 0. This implies that the corresponding 
projected wavefunction IV-") = V\ipo) satisfies S^\ip) = 
and {'4}\Si.\4>) ={]. 

As an aside, it is interesting to express Hq directly in 
terms of /^-fermions; it takes the form 



(r,r') 



(27) 



This Hamiltonian combines an imaginary spin-dependent 
hopping, with an imaginary pairing of the up-spin 
fermions only. 

Returning to the main task at hand, we define Majo- 
rana fermions as follows: 



i=V = - [iSr 



+ (T tr 



or, equivalently. 



= fri + fri 
^r — frf ~^ /rf 



(28) 

(29) 
(30) 
(31) 
(32) 



We note that the same mapping was recently employed in 
Ref. 34 to study Kitaev's honeycomb lattice model using 
5=1/2 partons. These objects satisfy the anticommu- 
tation relations 



{sr.,4,} = 0. 



(33) 
(34) 
(35) 



Moreover, from the form Eq. (28) it is clear that Sr is 
a singlet under projective spin rotations, while tr trans- 
forms as a vector. 

Expressing Hq in terms of Majorana fermions we have 



He 



-El 



iXrr'SrSr' + iXrr'^r ' U 



{r,r') 



where xt- 



Xrr' 



Xlr'/^ + 3x'w/2 and 



Xt'T"' 



(36) 



72- 



,/2. From this form is it clear that we simply 
have a theory of decoupled singlet and triplet Majorana 
fermions. Moreover, as long as x'^-r' 0) the s-fermions 
and t-fermions have different spectra. 
We also have the expressions 



1 



5; = --{iSrf. 



We note that GJ, = if and only if 



(37) 
(38) 

(39) 

(40) 



This is a Z2 form of the gauge constraint because the op- 
erator tl.t^t^Sr has eigenvalues ±1, and is the constraint 
appearing in the solution of Kitaev's honeycomb lattice 
model. 

In Appendix B, we show that IGG = Z2. This is done 
by directly proving that the only SU(2) gauge transfor- 
mations leaving Hq invariant are 17^. = 1 for all r, and 
Ur = —I for all r. It is also shown in Appendix B that 
space group and time reversal operations leaving Hq in- 
variant must be of the form 



S : Fr 



F, 



S(r) 



T:Fr^ {i(j^)Fr{i(J^)T:i , 



(41) 
(42) 



where tt^ and -kJ take values ±1 as a function of lattice 
site r. 

Now, we shall write down the low-energy effective 
gauge theory describing a Majorana spin liquid state. We 
note that essentially the same construction has been used 
previously for other Z2 spin liquid states.^'' For simplic- 
ity of notation, we consider an ansatz with only nearest- 
neighbor bonds. On every nearest-neighbor bond we 
place an Ising degree of freedom with a two-dimensional 
Hilbert space, acted on by Pauli matrices a^^, = a^,^ 
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and a^^, = cr^r^. is the Z2 vector potential, and is 
the Zi electric field. The local constraint becomes 

tltltlsr= n ^-r', (43) 

r'n.n.r 

where the product is over sites r' that are nearest- 
neighbors of r'. The Hamiltonian is 

{rr') 

-hJ2^rr'-Kj2 n <r'- (44) 
{rr') P rr'ep 

Wc take h,K > 0. The sums in the first two terms are 
over nearest-neighbor bonds. In the last term, the sum 
is over lattice plaquettes labeled by p, and the product 
rirr'Gp ^ product ovcr bonds in the perimeter of p. 

It is highly nontrivial to obtain the ground state phase 
diagram of H; there may well be a variety of phases. 
However, the physics is simple when K is sufficiently 
large, which is where the spin liquid phase arises. In 
this limit the Z2 gauge field enters its deconfined phase, 
where fluctuations of the Z2 magnetic field Ylrr'ep'^rr' 
are suppressed. The deconfinement is a robust property 
associated with a gap to Z2 vortex excitations, which 
are plaquettes where Yirr'ep'^rr' ~ other im- 

portant quasiparticle excitations are the fermions them- 
selves, which carry the Z2 electric charge as evident from 
their minimal coupling to the Z2 gauge field in Eq. (44). 
As for any state with a deconfined Z2 gauge field, this 
state is characterized in part by its Z2 topological order, 
as discussed for example in Ref. 35. 

Another important limit of H arises when h dominates 
over the other parameters. To be concrete, we set K = 
and assume h » |x*Mx*l- First setting = X* = 0) 
there is an extensively degenerate manifold of ground 
states, consisting of all states satisfying a^^, = 1 and 
the local constraint tl.t^t^Sr = 1- This is precisely the 
Hilbert space of the microscopic spin model. The de- 
generacy can be resolved using standard degenerate per- 
turbation theory in x* and x*- The first non- vanishing 
contribution occurs at second order, resulting in the ef- 
fective Hamiltonian 

HeS = ^ Jrr'Sr ■ Sr' , (45) 
{rr') 

where 

J^\rr'^^rr' ~^Xrr') fA£i\ 
rr' = • (4b j 

It is interesting to note that Jrr' is only antiferromag- 
netic when Xrr'iXrr' + Xlr') > 0. This suggests that 
states with negative xlr' ixlr' + Xrr' ) < are not likely 
to occur in microscopic Hamiltonians dominated by an- 
tifcrromagnctic exchange (but might reasonably occur if 
multi-spin exchanges are dominant). 



Wc close this section with some brief comments on 
projected wavefunctions for Majorana spin liquids, which 
may be obtained by applying the usual Gutzwiller projec- 
tion operator P to the ground state \tpo) of Hq- To sim- 
plify the discussion, we consider a state where x" = Xrr' 
and X* Xrr' ^^'^ non-zero only for r and r' nearest 
neighbors. When x"* = X*, then IGG = SU(2), and the 
projected wavefunction l'^) = Vlipo) can be associated 
with a low-energy effective SU(2) gauge theory, keeping 
in mind the caveats mentioned at the end of Sec. II. 
Now, if X*/x* is changed continuously from unity, the 
wavefunction does not change at all, because the pre- 
projected ground state \tpo} does not change. This oc- 
curs because \ipo) is a product of s-fermion and f-fermion 
ground state wavefunctions. It is interesting to note that 
this is so even though the IGG of the mean-field state 
- and hence the gauge group of the low-energy effective 
theory - is now Z2, as long as x*/x* 7^ 1- This is a rather 
dramatic illustration of the problematic association be- 
tween low-energy effective gauge theories and projected 
wavefunctions (sec Sec. II). 

In order to obtain a distinct Majorana spin liquid pro- 
jected wavefunction, we need to change Hq in such a way 
that IV'o) becomes different from the x* = X* ground 
state. One possibihty is to change the sign of x^/x*) as 
IV'o) does change when X^/x* crosses through zero. This 
only gives a single new wavefunction, as |V'o) is the same 
for all x*/x* < 0- Another possibility is, rather than sim- 
ply varying the ratio x'/xS to add further- neighbor hop- 
ping for, say, the s-fermions and not for the t-fermions. It 
will be interesting to study such wavefunctions in future 
work. 



V. MAJORANA SPIN LIQUID PROJECTIVE 
SYMMETRY GROUPS 

Here, wc show that, on any lattice, there is a one-to- 
one mapping between Z2 PSGs for Majorana spin liq- 
uids (Majorana PSGs), and SU(2) PSGs. Since classi- 
fication of SU(2) PSGs has already been done for some 
lattices, we can exploit those results to give a classifi- 
cation of Majorana PSGs. In this section, we first es- 
tablish the mapping between Majorana PSGs and SU(2) 
PSGs (Sec. V A). Next, we enumerate the four Majorana 
PSGs on the square lattice (Sec. VB). Finally, we con- 
sider frustrated mean- field ansatze of the form Eq. (36). 
We say an ansatz is frustrated if it contains at least one 
closed loop with an odd number of bonds, so that xt-r' 
(or Xrr') is nonzero for each bond in the loop. For in- 
stance, ansatze with triangular plaquettes are frustrated. 
In Sec. V G, we explain that any frustrated ansatz breaks 
time-reversal symmetry, and also note that there are no 
time-revcrsal-symmetric Majorana PSGs on the isotropic 
triangular lattice. 
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A. Correspondence with SU(2) projective 
symmetry groups 

Wc now show that, on any lattice, there is a one-to- 
one correspondence between Majorana PSGs and SU(2) 
PSGs. Precisely, we prove the following two statements: 
(1) Given a Majorana PSG, there is a corresponding 
unique SU(2) PSG, where the symmetry operations are 
realized exactly as in the Majorana PSG. (2) Given a 
SU(2) PSG, one can transform to a gauge where the sym- 
metry operations are realized exactly as they are in a 
unique corresponding Majorana PSG. In fact, the gauge 
needed in (2) is precisely the gauge used to classify SU(2) 
PSGs in Ref. 19. 

Wc begin by showing (1). A Majorana PSG is com- 
pletely specified by the following collection of transfor- 
mations under space group operations S, time reversal 
T, spin rotations, and Z2 IGG operations: 

S :Fr^Tr^Fs(r) (47) 

T:Fr^ TTjii(T^)Fr.{ia^) (48) 

Spin rotation : Fr UF^U'^ (49) 

IGG : Fr ±Fr. (50) 

Here, 7r^,7r^ = ±1 take values ±1 as a function of r. 
There is also an ansatz of the form Eq. (26) invariant 
under these transformations. It is shown in Appendix B 
that the forms Eq. (47) and Eq. (48) are the most general 
forms possible for an ansatz of the form Eq. (26). Such a 
set of transformation laws can be mapped into an equiv- 
alent set {i.e. same PSG) under a gauge transformations 
of the Z2 form Fr WrFr, where iTr takes values ±1 as 
a function of r. 

Now, we are free to continuously change Xrr' ^^'^ 
Xrr'i long as we change them in the same fashion on 
symmetry- related bonds. In particular wc can continu- 
ously tune Xrr' ^^^^ bonds, thus obtaining an 
ansatz with SU(2) IGG. This ansatz is clearly invariant 
under the same set of PSG transformations given above, 
and we have thus shown (1). 

Next, we show statement (2). Suppose we have a SU(2) 
PSG [with spin rotation symmetry realized naively ac- 
cording to Eq. (20)], where the IGG is generated by 

G' = ^Ml^Gi. (51) 

r 

Since G* must obey a SU(2) Lie algebra, wc know from 
Appendix A that either Mr G S0(3) or Mr = 0. We 
can make a gauge transformation so that Mr is either 
the identity matrix or zero. In fact, none of the Mr can 
be zero. Suppose Mr = for some site r, and is nonzero 
for some other site r'. In this case, there is no fermion 
bilinear joining r to r' that is invariant under the IGG. 
Since we want to consider only fully connected ansatze, 
we must then have Mr = 1 on all sites. (M = on 
all sites is also possible, but this would mean the IGG 
is not SU(2).) Since Mr S S0(3) for all r, the gauge 



transformation needed to turn into G* is unique, up 
to multiplication by an arbitrary gauge transformation in 
the Z2 center of SU(2). Upon completing the mapping 
to a Majorana PSG, this Z2 freedom will correspond to 
the Z2 gauge freedom to map one set of Majorana PSG 
transformations into another equivalent set. 

We have now gone to a gauge where [i?o,G*] = 0, 
where G' = 'l2r^r- By assumption, we also have 
[Hq.S"^] = 0. The most general quadratic Hamiltonian 
with these symmetries is 

Ho=tJ2 Xrr'tv{FrFl,). (52) 

{r,r') 

Next, consider a space group operation which acts 

by 

S:Fr^ Fsir)U^. (53) 

Because wc can multiply this transformation by any el- 
ement of the IGG, we are free to choose Urf^ = 1 for 
some arbitrary site vq. Now consider a site rg, joined 
to To by the bond {rQ,r'^). The gauge transformation 
Ur must transform the Hamiltonian on this bond into 
another bond Hamiltonian of the same form [as given in 
Eq. (52)]. This only happens if U^, = ±1. This conclu- 
sion holds for the whole lattice, because by assumption 
we can connect any site r to ro by some path of nonzero 
bonds. Therefore, we have 

S:Fr^w^Fs^r), (54) 

where tt^ takes values ±1 as a function of r. 

We proceed in essentially the same fashion for time 
reversal, which acts by 

T:Fr^ {ia^)FrUj. (55) 

For an arbitrary site ro, we choose Uj^ = (ia^). By the 
same argument as above, on all other sites we then have 
Uj^ = ±{ia^), leading to the desired result 

T:Fr-^ TrJii(T^)Fr{i(7^), (56) 

where irj takes values ±1 as a function of r. 

Finally, we note that by making appropriate IGG 

transformations, wc can choose spin rotations to act in 
the projective form iv UFrW. We have thus trans- 
formed to a gauge where the symmetry operations specify 
a unique Majorana PSG, and have shown statement (2). 

B. Majorana projective symmetry groups on the 
square lattice 

The square lattice space group is generated by the fol- 



lowing operations: 

)^{r^ + l,ry) (57) 

Ty:{r.,,ry)^{r,,ry + 1) (58) 

Px ■ {rx,ry) {-r^,ry) (59) 

Pxy ■■ {rx,ry) ^ {ry,r^). (60) 
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To specify a Majorana PSG, it is enough to specify the 
action of these operations, as well as time reversal, on 
the fermion operators. 

Exploiting the mapping between SU(2) and Majo- 
rana PSGs, and exploiting the results on classification 
of SU(2) PSGs on the square lattice in Ref. 19, we find 
there are four Majorana PSGs on the square lattice. We 
call these MAI, MA2, MBl and MB2. Ref. 19 refers to 
the corresponding SU(2) PSGs as SU(2)AnO, SU(2)A0n, 
SU(2)BnO, SU(2)B0n, respectively. To specify each of 
these PSGs, it is enough to give the action of the symme- 
try action on the singlet Sr fermions. The triplet tj, obey 
identical transformation laws. In Sec. VI, we give exam- 
ple mean-field states obeying each of these four PSGs. 



The MAI PSG is specified by: 

Ty . -^rx.ry ^ -^Tx ,ry-\-l (^^) 

Px ■ ^VxiTy ^ ( 1) rx,r-y (^^) 

T:s.,,.,^(-l)('-^+'^''^s,,,.,. (65) 
The MA2 PSG is specified by: 

Tx ■■ Sr^^ry Sr^ + l,ry (66) 

-^y • ^rx,ry ~^ ^rx,ry + l (67) 

Px ■■ Srx,ry {-lYyS-rx,ry (68) 

Pxy • '^r^.r-y ^ ^ry.r^ (^^) 

r:Sr^^,r^^i-l)^^-^+^y^Srx,ry. (70) 

The MBl PSG is specified by: 

Tx : Srx,ry ~^ ^r^+l,ry (71) 

Ty : Srx,ry ^ (" 1)"^ S,,,., + l (72) 

Px '■ ^rx,Ty ^ ( 1) '^^—Tx^Ty (73) 

Pxy • ^rx,ry ~^ ("1) ^^ry,rx (74) 

T:Srx,ry^{-l)^''^+'--^Srx,ry. (75) 

The MB2 PSG is specified by: 

Tx : Srx,ry ~^ ^rx + l,ry (76) 

Ty ■■ Srx,ry ^ {-lY" Srx,ry+1 (77) 
P:r ■■ .Sr.,r„ ^ (- 1)''« .S_., (78) 

Pxy ■ Srx,ry ^ (^1) "■^r^.r^ (79) 

T:Srx,ry^{-l)^'^-+'-^^Srx,ry. (80) 



C. Frustrated ansatze and time- reversal symmetry 
breaking 

It is worth noting that, for all the square lattice Majo- 
rana PSGs, the form of time reversal constrains the hop- 
ping to have a bipartite structure. That is, Xrr' ^^'^ Xrr' 
can only be nonzero if one of r, r' lies in the A sublattice, 
and the other site lies in the B sublattice. Therefore, on 
the square lattice, any frustrated mean-field ansatz of the 



form Eq. (36) breaks time reversal symmetry. [We say 
an ansatz is frustrated if it contains at least one closed 
loop with an odd number of bonds, so that Xrr' Xrr') 
is nonzero for each edge in the loop.] 

In fact, this statement holds on any lattice for Majo- 
rana spin liquid ansatze. Suppose we have a frustrated 
ansatz where Xrr' 7^ on the bonds of an odd-length 
closed loop. (We could just as well consider a loop with 
xtr' 7^ 0-) The time reversal operation takes the general 
form given in Eq. (42); on the Sr fermions we have 

r : St- -5- Trjsr- (81) 

If we choose ttJ = 1, this transformation simply sends 
Xt-r' ~^ ~Xrr'- ^or loops of cvcn length this change can 
be compensated by an appropriate choice of gauge trans- 
formation n^, but for a loop of odd length this is im- 
possible, so any frustrated ansatz breaks time reversal 
symmetry. 

To illustrate the strong restriction this imposes on Ma- 
jorana spin liquid ansatze, we consider the isotropic tri- 
angular lattice with 6-fold rotation symmetry and full 
translation symmetry of the triangular lattice. Using 
only these symmetries, it is easy to show that any ansatz 
contains a closed loop of length three (this is true even 
if the ansatz has vanishing nearest-neighbor hopping). 
Therefore, if one insists on time reversal invariance, there 
are no Majorana PSGs on the isotropic triangular lattice. 

VI. MAJORANA SPIN LIQUIDS ON THE 
SQUARE LATTICE 

In this section we study the simplest ansatz for each 
of the four Majorana PSGs on the square lattice. The 
ansatze we consider are the simplest in the sense that 
they include only the shortest-distance hopping permit- 
ted by symmetry. In the case of MAl and MBl PSGs, 
this is nearest-neighbor hopping, while for MA2 and MB2 
PSGs, we have fourth-neighbor hopping. Each state is 
associated with a corresponding mean-field SU(2) spin 
liquid, obtained by setting x'^ — X*- For the MAI state 
this is the uniform resonating valence bond (RVB) state, 
and for the MBl state this is the 7r-flux state. 

For each state, we study the mean-field excitation spec- 
trum. In all cases except the MBl state, we find a nested 
Fermi surface, which we expect to give rise to instabili- 
ties upon going beyond mean-field theory and including 
short-range interactions of fermions. However, at least 
for the MAI and MB2 states, these instabilities are ex- 
pected to be logarithmic in nature (similar to BCS in- 
stability), and therefore may play a role only at very low 
temperatures. The MA2 state has discrete points with 
a 2 = 4 excitation spectrum {i.e. energy goes like the 
fourth power of momentum) , which arc likely to lead to a 
stronger instability that will be important at higher tem- 
peratures. The likely consequence of these instabilities is 
magnetic order, but it should be noted that the resulting 
ordered states are exotic, supporting gapped Z2 vortex 
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excitations and gapped fcrmionic spinons. This occurs 
because deconfinment of the gauge field is protected 
by the Z2 vortex gap, and is thus robust to arbitrary 
small perturbations. 

In contrast to the other states, the MBl state has gap- 
less Dirac points, and is thus dubbed the MBl-Dirac 
state. The MBl-Dirac state is stable to arbitrary small 
perturbations, provided that space group, time rever- 
sal, and spin rotation symmetries are respected. This 
is shown in Sec. VII, where some other properties of the 
MBl-Dirac state are also discussed. 



MAI state 



The simplest mean-field ansatz with MAI PSG has the 
Hamiltonian 



r r 

When X* = x*j this Hamiltonian reduces to pure imagi- 
nary hopping of the S = 1/2 fra fermions, with the hop- 
ping phases such that there is zero magnetic flux through 
each plaquette. Therefore, we obtain the SU(2) uniform 
RVB state at this special point. 

Since both s- and t-fermions have the same spectrum, 
only with different coefficients, it is enough to focus on 
Hqs, the s-fermion part of Hq. Diagonalizing Hqs, we 
find the single-particle spectrum is given by 



Es{k) = 4|xs|| sinfc^ -|- sinfcyl, 



(83) 



where — tt < kx,ky < tt, with the proviso that k and —k 
are equivalent points due to the Majorana nature of the 
fermions. 

There are lines of gaplcss excitations for ky = —k^ and 
ky = ±7r -|- kx, which constitute the familiar diamond- 
shaped Fermi surface of the uniform RVB state. More- 
over, it is straightforward to show that invariance under 
MAI PSG transformations requires these lines to be gap- 
less at the quadratic level. This means, without breaking 
some symmetry, it is impossible to remove the Fermi sur- 
face nesting by adding further neighbor hopping. Upon 
going beyond mean-field theory and incorporating short- 
range fermion interactions, it is natural to expect that the 
nested nature of the Fermi surface leads to an instabil- 
ity to Neel magnetic order. This is well-known to occur 
when x" = X*; ^-iid we have verified it more generally, 
treating an on-site interaction in mean-field theory, and 
finding an instability to Neel order for arbitrarily small 
interaction.'^^ 



B. MA2 state 

The simplest mean-field ansatz with MA2 PSG has the 
following singlet part of the Hamiltonian: 



Os 



Hqs 



^x. 



n 



—x+2y SrSr—2x+y 



(84) 



For this PSG, first, second and third neighbor hoppings 
are required to vanish, but fourth-neighbor hopping is 
allowed by symmetry. It is important to remember that 
hopping of Majorana fermions carries an orientation; the 
pattern of orientations for the fourth-neighbor hopping of 
Hqs is illustrated in Fig. 1. The single-particle spectrum 
is given by 



Es{k) = 4|xs| sm{2kx + ky) + sm{kx + 2ky) 



— sva.{kx — 2ky) + sm{2kx — ky) 



(85) 



where — tt < kx,ky < tt (again remembering that k and 
—k are equivalent points). There arc; lines of gaplcss ex- 
citations for ky = —kx, ky = ±7r -|- kx, ky ~ ±7r/2, and 
kx = ±7r/2. It can be shown that these gapless lines are 
protected (at the quadratic level) by MA2 PSG trans- 
formations. Moreover, we again expect an instability to 
magnetic order due to Fermi surface nesting. 

It is interesting to note that the spectrum at A; = 
{n/2, —■ 7r/2) has a 2; = 4 character. Expanding to lowest 
order near this point, we have 



Es[k~{7T/2,-7T/2)]^8\x'\\kxkl 



k^ ky\ 



(86) 



While we have not studied the issue in detail, the strong 
infrared singularities from such a point may lead to a 
strong instability to magnetic order (stronger than the 
usual logarithmic instability arising from Fermi surface 
nesting) . We also note that we have not verified whether 
the z = 4: nature of this point is protected by MA2 PSG 
transformations. 



C. MBl state 

The simplest mean-field ansatz with MBl PSG has the 
following singlet part of the Hamiltonian: 



H, 



Os 



R 



SRiSR2 + SriSR+x,! 



+ Sii2SR+2y,l — Sii2Sr+x,2 



(87) 



We use a two-site unit cell as shown in Fig. 2, where 

the orientations of the nearest-neighbor hopping are also 
shown. Sites are labeled by pairs (H, i), as described 
in the caption of Fig. 2. When x* = X*) this Hamilto- 
nian reduces to pure imaginary hopping of the 5* = 1/2 
fra fermions, with the hopping phases such that there 
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FIG. 1. Orientation of fourth-neighbor hoppings Xrr' in the 
MA2 state. All other hopping amplitudes can be obtained by 
a translation of those shown. 
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FIG. 2. Orientation of nearest-neighbor hoppings in the MBl 
state. The two-site unit cell is indicated by the dashed line. 
Unit cells are labeled by ii = UxX + 2nyy. Sites are labeled 
by the pair (R, i), where i = 1 corresponds to the lower site, 
and the upper site is i = 2. 



is magnetic flux of tt through each plaquette. Therefore, 
we obtain the SU(2) tt-Aux state at this special point. 

The single-particle spectrum is two-fold degenerate, 
and is given by 



Es{k) = 2|x1y^4-2cos(2A;^) -2cos(2fcy), (88) 



where — TT < fcj; < TT, — 7r/2 < ky < tt/2, and again 
it should be remembered that k and —k are equivalent 
points. There are gapless Dirac nodes at fc = (0, 0) 
and k = (7r,0), so we dub this state the MBl-Dirac 
state. A low-energy theory for this state is derived in 
Appendix C, and physical properties of the state are dis- 
cussed in Sec. Vll. In particular, it is shown that the 
MBl-Dirac state is a stable spin liquid phase. 



FIG. 3. Orientation of fourth-neighbor hoppings in the MB2 
state. The two-site unit cell is indicated by the dashed line. 
Unit cells are labeled hy R = UxX + 2nyy. Sites are labeled 
by the pair {R,i), where i = 1 corresponds to the lower site, 
and the upper site is i = 2. 



D. MB2 state 

The simplest mean-field ansatz with MB2 PSG has the 
following singlet part of the Hamiltonian: 



R 



— Sr\Sr-2x,2 + SR2SR+2x+2y,l — SR2SR+x+2y,2 



SR2SR-x+2y,2 — S R2S R-2x+2y ,1 



(89) 



We again use a two-site unit cell as shown in Fig. 3, where 
the orientations of the fourth-neighbor hopping arc also 
shown. As in the MA2 state, first, second and third 
neighbor hopping is forbidden by symmetry, but fourth- 
neighbor hopping is allowed. 

The single-particle spectrum is two-fold degenerate, 
and is given by 

Es{k) ^ 8V2\xs\\cos{kx)cos{ky)\ 

X y2^^cos(2fc^y^^cos(2fc^, (90) 



where — tt < kx < tt, — 7r/2 < ky < 7r/2, and again 
it should be remembered that k and k are equivalent 
points. There are gapless Dirac nodes at fc = (0, 0) and 
fc = (tt, 0), as well as gapless lines for kx = ±7r/2. It can 
be shown that the gapless points and lines are protected 
(at the quadratic level) by MB2 PSG transformations. 
We again expect an instability to magnetic order due to 
Fermi surface nesting. 



VII. PROPERTIES OF MBl-DIRAC STATE 

Here, we briefly discuss the low-energy effective theory 
of the MBl-Dirac state, and some of its physical proper- 
ties. In particular, we show that the MBl-Dirac state is a 
stable phase, and discuss the effect of Zeeman magnetic 
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field. We also discuss some nearby gapped symmetry- 
breaking phases, and show that states where vortices 
are bound to an odd number of Majorana fermions (and 
thus have non-Abelian statistics) can occur when VT- 
breaking order coexists with columnar dimer order. 

The continuum low-energy theory (at the mean-field 
level) is worked out in Appendix C. This theory is ob- 
tained by linearizing about the two gapless Dirac nodes 
at fc = (0,0) and h = (tt, 0). The imaginary-time action 
S = J dT(Pr£ is specified by the Lagrangian density 

c = ni^^d;]^ + ¥[i^,d'^W. (91) 

Here, 4' is a four-component real fermion field arising 
from the singlet s-fermion on the lattice. For each triplet 
fermion V, is the corresponding four-component con- 
tinuum field. We denote the 2x2 Pauli matrices acting 
in the two-dimensional space of each node by r*, and 
/i* Pauli matrices act in the flavor space mixing the two 
nodes (see Appendix C for more detail). The space-time 
index /U = 0, 1,2, and we define 7^ = (t^,t^,— r^). The 
singlet (triplet) fermions have velocity Vg (vt), which en- 
ter via the derivatives = {do,Vsdi,Vsd2), and 9* = 
{do,Vtd\,Vtd2)- Finally, we define ^ = ^^(— ir^), and 
similarly for 

This low-energy theory is in fact stable to the addition 
of small perturbations beyond mean-field theory, and the 
MBl-Dirac state is thus a stable phase. Because excita- 
tions of the deconfined Z2 gauge field are gapped, cou- 
pling to the Z2 gauge ficild has no effect on the fermions 
at low energies and can be safely ignored. While the re- 
sulting low-energy theory does not capture the Z2 topo- 
logical order that is present, it does correctly describe 
the universal behavior of correlation functions of local 
observables. Moreover, it should be noted that ignoring 
fermion-gauge field coupling means that our low-energy 
theory if we add perturbations large enough to desta- 
bilize the MBl-Dirac state is only capable of describing 
phases where the Z2 gauge field remains deconfined. 

We also need to consider perturbations involving only 
the fermion fields. The action S is invariant under the 
renormalization group (RG) scale transformation r — ^ 
e^T, r e^r, ^ e~^^, and e~^$'. This 

amounts to the statement that we have an RG fixed point 
with dynamic critical exponent z = 1, where the fermion 
fields have unit scaling dimension. Perturbations that are 
relevant or marginal under this RG transformation can 
destabilize the phase; such perturbations are fermion bi- 
linears with no derivatives (dimension two — relevant), 
and with one derivative (dimension three - marginal). 
Making use of the symmetry transformations given in 
Appendix C, it can be shown that all bilinears with no 
derivatives are forbidden by the combination of spin ro- 
tation, space group and time reversal symmetries. The 
only single-derivative bilinears allowed by symmetry are 
shifts of the velocities Vg and Vt- All other perturbations, 
including quartic interactions of fermions, are irrelevant 
under the RG. Therefore the MBl-Dirac state is a stable 
phase. 



h=0, y/Ar'=0.5 
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FIG. 4. Plot of the h = Q singlc-particlc excitation spectrum 
for the MBl-Dirac state, along the line ky = 0, from kx = Q 
to kx = TT, for x*/x* = 0-5- The solid line shows the six-fold 
degenerate t-fermion energy, while the dashed line is the two- 
fold degenerate s-fermion energy. The gapless Dirac nodes 
are evident at fcx = and fcx = tt. 



h//=0.8, y/Ar'=0.5 
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FIG. 5. Single-particle excitation spectrum of the MBl-Dirac 
state, for /i/x* = 0.8 and X^/x* = 0.5, along the hnc ky = 0. 
Here, the solid lines are the energies in the t^-t^ sector, and 
the dashed lines the energies in the s-t^ sector. Around the 
nodes, small Fermi pockets are opened in both sectors. Away 
from the nodes, the t-fermion spectrum is split into three 
distinct branches. 



The same scaling considerations described above also 
imply the heat capacity C{T) oc and the magnetic sus- 
ceptibility x{T) T. Correlation functions of fermion bi- 
linears fall off as 1/r^ in space and 1/t** in time. This im- 
plies in particular that k = (0, 0), (tt, 0), (0, tt) and (tt, tt) 
spin correlations fall off with these power laws. 

We now return to the lattice to discuss the effect of 
Zeeman magnetic field. In the low-energy effective the- 
ory, a z-axis Zeeman field gives the following contribution 
to the Hamiltonian: 

r r 

where Cst and Ctt are dimensionless, non-universal con- 
stants. Two terms appear in this Hamiltonian because 

each of these terms has precisely the same symmetries as 
the microscopic magnetic field operator . In the 
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h//=0.8, Ar'//=0.5 




FIG. 6. Detail of the single-particle excitation spectrum of 
the MBl-Dirac state, along the line ky = 0, near the node at 
kx = 0. The solid lines arc the energies in the t^-t^ sector, 
and the dashed lines the energies in the s-iP' sector. The 



parameters are the same as in Fig. 5, that is /i/x* 
XVX* = 0.5. 



h//=0.8, y/Ar'=-o.5 



0.8 and 




FIG. 7. Detail of the single-particle excitation spectrum of 
the MBl-Dirac state for x7x' = -0.5 and /i/x* = 0.8. The 
plot is again along the line ky = 0, near the node at kx = 0. 

The solid lines are the energies in the t^-t^ sector, and the 
dashed lines the energies in the s-t^ sector. It is evident that 
a gap has opened in the s-t^ sector. The spectrum away from 
the nodes is very similar to the case Xa/xt = 0.5, as plotted 
in Fig. 5. 



presence of one of these terms (say Cgt = but Cu 7^ 0) 

and short-range fermion interactions, we expect the other 
term will be generated. Therefore, both terms must be 
included. 

We focus first on the sector of the theory involving s- 
and t^-fermions. Neglecting coupling to the gauge field 
and fermion interactions (which, for the effect of a small 
Zeeman field on the low-energy spectrum, is not an ap- 
proximation), the single-particle excitation spectrum is 
given by 



e5!(fc) = ix' + x')h x'Yfl + h% , (93) 

where we have defined hst = hCst/4: and 



/fe = J4 - 2 cos(2A;^) - 2 cos{2ky). (94) 



Each energy level is two-fold degenerate, and it should 
be recalled that that k and k are equivalent points. 
The spectrum in the sector involving t^- and t^-fermions 
can be obtained from Eq. (93) by putting X* and 

hst htt = hctt/4:. 

The effect of the Zeeman field differs considerably de- 
pending on whether Xs and xt have the same or opposite 
signs. When Xs/xt > 0, a small Zeeman field opens 
Fermi pockets around the Dirac nodes in both the s-t^ 
and t^-t"^ sectors. On the other hand, when Xs/Xt < 0, a 
small Zeeman field opens a gap in the s-t^ sector. This 
is illustrated in Figures 4-7 (in all cases shown wc choose 
Cst = Ctt = !)• These spectral features will be manifest 
in the spin structure factor S{q,io), which can be mea- 
sured by neutron scattering. Calculation of S{q,uj) for 
the MBl-Dirac state is left for future work. We note that 
similar interesting behavior of S' = 1 Majorana spinous 
in a Zeeman field was also found in the exactly solvable 
model of Ref. 14. 

Finally, returning to the continuum effective theory, 
we consider the properties of some gapped phases nearby 
to the MBl-Dirac state. In particular, given the appear- 
ance of non-Abclian statistics in the _B-phasc of the Ki- 
taev honeycomb lattice model, as well as in one of the 
exactly solvable models with 5=1 Majorana spinous, 
it is natural to ask whether such statistics can also occur 
upon opening a gap in the MBl-Dirac state. We an- 
swer this question in the affirmative below, although - at 
least if we maintain spin rotation symmetry - it seems to 
be necessary both to induce a dimerization of the square 
lattice, and to break time-reversal symmetry. It is likely 
that there are other routes to non-Abelian statistics if 
we allow for breaking of spin rotation symmetry, but we 
have not studied this possibility. 

In the gapped phases we consider, various symmetries 
are broken - this symmetry breaking can be induced ex- 
plicitly, or could potentially be induced spontaneously as 
a result of sufficiently large fermion interactions. In all 
these phases, Z2 vortex excitations remain gapped and 
the Z2 gauge field remains in its deconfined phase. In 
particular, we consider the following perburbation to the 
Lagrangian: 

(95) 

For simplicity of discussion, we restrict attention to the 

case mfj , m,(j , m^jj , m]j > 0. As discussed in Ref. 16, the 
Z2 vortices in a system such as this one have non-Abelian 
statistics when the Chern munbcr z/ is odd, stemming 
from the binding of an odd number of Majorana modes 
to each vortex. 

First wc consider the case m^, 7^ and mfj, m^, = 
0, where the Lagrangian breaks parity and time reversal 
symmetries, while preserving the product of parity and 
time reversal, as well as other microscopic symmetries. In 
this case each species of fermion on the lattice has Chern 
number 1^ = 1, for a total Chern number v = A, giving rise 
to four gapless, co-propagating Majorana edge modes. 
(We note that if we allow ttiq and ttiq to have opposite 
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sign, we can also obtain v = 2.) Because the Chcrn 
number is even, the vortices do not have non-Abehan 
statistics. 

Next; wc consider m^, = and mfj, m^j, ^ 0. This 
corresponds to inducing a columnar dimcrization of the 
square lattice. Each lattice fermion has v = 0, so the 
total Chern number = 0, and the fermion spectrum in 
this state is topologically trivial. In particular there are 
neither gaplcss edge states nor non-Abelian statistics. 

Finally, we consider the case where all the mass terms 
are nonzero, which corresponds to co-existing columnar 
dimer and time-revcrsal-breaking orders. First we con- 
sider the s- fermion sector. As long as > ml^, then 

V = 1. When rriQ = mfj, the gap closes, and for 

< mfj wc obtain a topologically trivial state with 

V = Q. The same statements hold for the i-fermion sec- 
tor, except that the total i-fermion Chern number when 

> is V = i. Based on this, we note that if 
rriQ > but m\j < m^,, we have a total Chern num- 
ber v — I, giving rise to one gaplcss chiral Majorana 
edge mode. If we reverse the above inequalities, we in- 
stead have v = 3, with three Majorana edge modes. Non- 
Abelian statistics arise in both these cases. 



VIII. DISCUSSION 

In this paper, we showed that, in the fermionic parton 
approach to 5 = 1/2 spin liquids with SU(2) symmetry, 
spin rotations can be realized in two distinct ways. In 
addition to the familiar naive realization, where spin ro- 
tations act on the S' = 1/2 partons with no accompanying 
gauge transformation, we also found a projective realiza- 
tion, where SU(2) spin and gauge rotations are locked 
together. This projective realization leads to spin liquids 
with 5 = 1 and 5 = Majorana fermion excitations cou- 
pled to a deconfined Z2 gauge field. We discussed the 
projective symmetry group (PSG) classification of such 
states, showing that their PSGs are in one-to-one corre- 
spondence with SU(2) PSGs in the existing classification. 
To illustrate these results we studied states in each of the 
four Majorana PSGs on the two-dimensional square lat- 
tice, finding that one, the MBl-Dirac state, is a stable 
phase with gapless Fermi points. This phase exhibits in- 
teresting behavior in a Zeeman magnetic field, and has 
nearby gapped phases supporting non-Abelian statistics. 

We now discuss the challenging and very interesting 
question of where Majorana spin liquids may be found in 
models (beyond the exactly solvable models of Refs. 11, 
12, 14, and 15) and in real systems. In Sec. VC, we 
noted that frustrated mean-field ansatze for Majorana 
spin liquids break time reversal symmetry. It is likely 
that such time-reversal-breaking states will have lower 
energies in frustrated 5=1/2 magnets than Majorana 
spin liquids with time reversal symmetry, because they 
may be better able to gain exchange energy from all the 
bonds of the frustrated lattice. Because frustration is 
expected to be an important ingredient in stabilizing spin 



liquid phases, it may be more productive to search for 
time-revcrsal-breaking Majorana spin liquids than their 
fully symmetric cousins. It will be interesting to classify 
and study such states in future work. 

Another promising place to look for spin liquids in 
general is in weak Mott insulators; that is, just on the 
insulating side of a Mott metal-insulator transition. 
A given spin liquid state is naturally associated with a 
corresponding state, not necessarily unique, obtained by 
condensing charge-carrying excitations in the spin liquid. 
This state can be a metal, superconductor or band insu- 
lator. This may also be a promising place to look for 
Majorana spin liquids; indeed, a Majorana spin liquid 
with projective spin rotation symmetry has already been 
discussed in a continuum effective theory of the honey- 
comb lattice Hubbard model (phase B2 of Rcf. 37). 

Majorana spin liquids could be studied in Hubbard 
models using the SU(2) slave-rotor approach,^^"^" which 
is an extension to the Hubbard model of the SU(2) gauge 
theory of the Heiscnberg model. For single band H\ib- 
bard models at half- filling on bipartite lattices, in ad- 
dition to the SU(2) spin rotation symmetry, there is a 
SU(2) pseudospin symmetry, of which U(l) charge rota- 
tions is a subgroup. **^''*^ In the slave rotor framework, we 
expect that the state obtained from a Majorana spin liq- 
uid by condensing charge-carrying excitations will have a 
spontaneous locking of spin and pseudospin symmetries; 
that is, SU(2) x SU(2) symmetry is broken to the diago- 
nal SU(2) subgroup of simultaneous spin and pseudospin 
rotations. This expectation is corroborated by the re- 
sults of Ref. 37. (There, phase B2 is adjacent to phase 
B, which has such a SU(2) x SU(2) -J> SU(2) symmetry 
breaking.) In more realistic situations, the SU(2) pseu- 
dospin is broken down to U(l) charge rotations, and we 
may expect instead a locking of charge rotations to a 
U(l) subgroup of spin rotations, corresponding to break- 
ing SU(2) X U(l) ^U(l). 

On a different note, it is interesting to ask whether our 
starting point of 5 = 1/2 fermionic partons is necessary 
to describe the Majorana spin liquids discussed here - 
it is not. We could have introduced Majorana partons 
from the beginning, using the representation of the spin 
operator in Eq. (37) and the constraint Eq. (40). Con- 
sidering mean-field Hamiltonians of the form Eq. (36), 
we would again conclude that [as defined in Eq. (39)] 
and not 5' commutes with the Hamiltonian. The status 
of spin rotations in such a state is potentially a confus- 
ing issue, but the results of this paper - and the starting 
point of 5 = 1/2 partons - clarify that this state is in 
fact spin-rotation invariant. 

Actually there is some freedom in the expression of 5* 
in terms of Majorana partons; more generally, we may 
have 

1 i 

5* = -- [(1 - x)iSrtri + (1 + x)-eijktrjtrk] ■ (96) 

Setting X = gives the representation wc obtained from 
5=1/2 fermions, and x = — 1 is the representation 
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used by Kitaev to solve his honeycomb lattice model. -'^^ 
Putting a; = 1 gives S^. = T^, as in the representation 
developed in Refs. 43-45, and employed by Biswas et. 
al}^ (see Appendix D). Unlike in that representation, 
however, even when x = 1 there is still an s-fermion that 
enters via the constraint. 

We close with a discussion of some open issues for 
future study. It would be interesting if the formal- 
ism developed here can reproduce the S — I Majo- 
rana spin liquid phases of the exactly solvable models of 
Refs. 11, 12, 14, and 15. We have made some attempts in 
this direction for the model of Ref. 12, but so far have not 
been successful. It would also be interesting to study pro- 
jected wavefunctions for Majorana spin liquids, which we 
discussed briefly in Sec. IV. Finally, there are no doubt 
other circumstances where a continuous global symme- 
try is realized projectively in a spin liquid or other exotic 
phase. For example, suppose we consider a 5 = 1/2 sys- 
tem with only global U(l) symmetry, which we can think 
of as a system of strongly correlated bosons. Treating this 
system using S = 1/2 fermionic partons, the U(l) sym- 
metry can be realized projectively by a locking to a U(l) 
subgroup of the SU(2) gauge group. Due to the lower 
symmetry, we expect that the form of the mean-field 
Hamiltonian is less constrained, and that a wider vari- 
ety of spin liquid states may occur. Moreover, projective 
continuous symmetry can be realized in this system with 
bosonic S = 1/2 partons, where the gauge group is U(l), 
by a locking of gauge and global U(l) symmetries. These 
and other similar states may have unexpected properties, 
and may prove important in broadening the understand- 
ing of exotic states of matter. 

While this paper was being finalized, we learned that 
T. Senthil has independently obtained some of the main 
results presented here."^^ 
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Appendix A: Restrictions on Mr 

In Sec. Ill, it is asserted that, starting from the form 

T; = S"; + MiPGi (with M;^' an arbitrary real 3x3 
matrix), the requirement [T^, T^] = ie^^'^T^ implies either 
Mr = or Mr e SO (3). In this Appendix we prove this 
assertion. 

The required commutation relation implies the follow- 
ing ciquation for the matrix M (here and below, we drop 
the site label r): 

^ijkj^km ^ ]^ilj^jk^lkm_ ^^^^ 

We will now find all possible solutions of this equation 
for M. We proceed by making the singular value decom- 



position M = allDV, where a = ±1, U,V € S0(3), and 
D = diag(di, (i2, rfs), where di > 0. Making use of the 
identities WW'^ = 5^^ and U'"^' U"' U^^' e^'^'^' = e'J'=, 
and similarly for V, Eq. (Al) can be brought to the form 

^^ijm-Qmk _ jjin £)jm ^nmk (A2) 

This gives the equation adz = did2 and its two cyclic per- 
mutations. Clearly one solution is di = 0, corresponding 
to M = 0. If any one of the di is nonzero, then clearly 
they must all be nonzero. Moreover, in the case of a 
non-zero solution, since the di are positive, we must have 
a = 1. It is then trivial to show that the only non- 
zero solution is di = ^2 = ^3 = 1, which corresponds to 
M G S0(3). 

Appendix B: Invariant gauge group and form of 
symmetry operations 

Here, we show that Majorana spin liquid mean-field 
ansatze, with Hq as given in Eq. (26), have IGG = Z2. 
This statement holds as long as Xrr' on some bonds 
- we recall that we always assume this to be the case, 
since if Xrr' ~ everywhere there is no distinction be- 
tween naive and projective spin rotations, and we simply 
have an ansatz with IGG = SU(2). We also show that 
the space group and time reversal operations leaving iJo 
invariant take a simple form. 

To show that IGG = Z2, consider first any two sites r 
and r' for which x^^, is nonzero. The Hamiltonian for 
this bond is 

^rr' = ixlr' ^^^4') + ^xlr' tr(a*F,(T' fJ, ) • (Bl) 

We wish to find all gauge transformations Fr — 
FrUr leaving Hrr' invariant. It is useful at 
this point to note that the 16 fermion bilin- 

ears {tr(j;i^J,), tr(a*J^,Fj,), tr(i^,c7*Fj,), tr(a^F^a^Fj,)} 
form a complete, linearly independent set of all fermion 
bilinears connecting the two sites r and r'. In 
particular, this means that any gauge transformation 
leaving Hrr' invariant must separately leave the 
and terms invariant: gauge transformations send 
the x^ term to a linear combination of the bilinears 
{tr(FrF^,),tr(Fr(T'F^,)}, while the term becomes a 
linear combination in the distinct subspace of bilinears 
spanned by {tr(a*F^Fj,), tr(a*F^o-^Fj,)}. 

In order for the x^ term to be invariant, we must have 

Ura'Ul,=a\, (B2) 

for all i = 1,2,3. The only solution to this equation is 
Ur = = ±1. The x^ term is also invariant under such 
a gauge transformation, so we do not need to consider it 
separately. 

Next we need to go beyond the original two sites r 

and r' . Suppose there is a third site r" connected to r 
such that xir" 7^ or Xrr" 7^ 0- In order to leave Hrr" 



16 



invariant, It is easy to see that we must have Ur" = 
Ur- Because we assume the ansatz to be connected, we 
can repeat this procedure to determine Ur for the entire 
lattice. Therefore we have shown that the only gauge 
transformations leaving Hq invariant are Ur = 1 for all 
r, and Ur = —I for all r, and thus IGG = Z2. 

Now we discuss some restrictions on the form of sym- 
metry operations. As usual, we suppose that Hq is in- 
variant under a set of space group symmetries labeled by 
S, and also under time reversal symmetry T. As we see 
below, these operations must take the form 



S : Fr ^ Fs(r)'^r 

T:Fr^ {ia'^)Friia^)TTj, 



(B3) 
(B4) 



where tt^ and nj take values ±1 as a function of lattice 
site r. 

We first consider the space group transformations. The 
most general action of a space group transformation S on 
fermion operators is 



S : Fr Fsir)U; 



(B5) 



We can think of this operation as a composition of 
the gauge transformation Ur followed by the operation 
Fr Ps(r) ■ It is useful to note that the gauge transfor- 
mation Ur need not leave Hq invariant, but must trans- 
form the Hamiltonian on each bond Hrr' into another 
Hamiltonian of the same form, meaning that the result of 
the gauge transformation can be absorbed into a change 
of xlr' ^^'^ xlr'- Foi' two sites r and r' with xlr' 0; 
the only such transformations are U^ = U^, = ±1, and 
Ur = —U^i — ±1. Now, as above in determining the 
IGG, suppose a third site r" is connected to r such that 
Xrr" 7^ or Xrr" 7^ 0- Order for Hrr" to trans- 
form into another bond Hamiltonian of the same form, 
we must have U^,, = ±1. Following this procedure to 
extend the transformation to the whole lattice, we find 
Ur can only take values ±1 for all r, and thus Eq. (B3) 
holds. 

Next, we consider time reversal. The most general 
realization of time reversal is 



T:Fr^ {i(j'')FrUj {ia 



(B6) 



where the presence of (icr^) on the right is a convention, 
and it should be kept in mind that the transformation 
is anti-unitary. This transformation must leave Hrr' in- 
variant. Acting on Hrr', we have 

r : Hrr' ^ -ixlr' tr{F^Uj {Uj,)^ F^ 



rr' 



r,tr{F,Uj{Uj,)^Fl,) 
tv{a'FrUja\Uj,)^Fl,) (B7) 



As long as Xrr' 7^ 0, the only transformations leaving 
Hrr' invariant are Uj = —Uj, = ±1. We can follow 
the above procedure to extend Uj to the whole lattice, 
and find that Uj can only take values ±1 for all r, thus 
showing Eq. (B4). 



Appendix C: Continuum field theory for MBl-Dirac 

state 

Here, we work out the low-energy effective theory for 
the MBl-Dirac state, at the mean- field level, and quote 
the action of the microscopic symmetry operations in the 
low-energy fermion fields. As discussed in Sec. VII, this 
low-energy theory is also valid beyond mean-field theory. 

We begin by working with the s-fermion part of the 
mean-field Hamiltonian. Because the development for 
the f-fermions exactly parallels the treatment below, 
there is no need to go through it explicitly. We have 



X'^nT'ij{k)skiS-kj, 



(CI) 



where the sum is over the Brillouin zone Ifc^l < tt, \ky\ < 
7r/2, and 



m{k) 



(C2) 



The nodes are at fc = and k = K = (tt. 0). Letting q 
be small on the scale of the zone size, we have 



m{q) = 2q^T^ + 2gyr^ 
m{q + K) = -2q^T^ + 2gyT\ 

We define continuum fields by writing 



V'2 ~ T 



Sql 
Sq2 

1 f Sq+K,l 



(C3) 
(C4) 



(C5) 
(C6) 



where r' are 2x2 Pauli matrices. These fields obey the 
continuum Hamiltonian density 



(C7) 



where we have introduced the nodal velocity Vg ■ Defining 
V'l = i^ii^i'^^) S'nd 'tp2 = V'2'(— «T^), the imaginary-time 
Lagrangian density is 

Cs = M^lf^d;,]^! + Milf.df,]ij2, (C8) 
where /U = 0, 1, 2, 

d;, = {do,v,duv,d2), (C9) 

and 

j^ = {r',r\-T'). (CIO) 
We define the four-component spinor 



^1 

i>2 



(Cll) 
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and define /i* Pauli matrices acting in the flavor space 
mixing tpi and ^2- For example, 




(C12) 



For each triplet fermion t*, we proceed identically and 
introduce a four-component spinor field The full con- 
tinuum Lagrangian density is then 

c = ^i^^d;]^ + nh^.diw, (C13) 

where 

dl = {do,Vtd„Vtd2), (C14) 

and we have introduced the nodal velocity Vt of the t- 
fermions. 

The action of space group and time reversal symme- 
tries on the lattice fermion fields is given in Eqs. (71-75). 
It is straightforward but somewhat tedious to work out 
the action of these operations on the continuum fields. 
As these symmetries act identically on 4* and we now 



quote their action on ^: 

: *(r) ^ /z3*(r) (C15) 

Ty : ^ /i^*(T-) (C16) 

: *(r) ^ TV'*(r') (C17) 

P^cy : *(r) -^l{r'+ t'){h' + /x^)*(r') (C18) 

r : *(r) (ir2)(V)*(r). (C19) 



Here, r = {rx,ry), : r ^ r' = {-rx,ry), and P^y : 
r ^r' = {ry,rx). 

Appendix D: Approach of Biswas et. al. 

A different approach to S = 1 Majorana spin liquids 
was introduced by Biswas, Fu, Laumann and Sachdev 
(BFLS),^'"* where a mean-field theory of a particular Ma- 
jorana spin liquid on the triangular lattice (BFLS state) 
was constructed. This approach has some similarities 
to ours, but differs in the absence of a singlet s-fermion. 
Because we feel it may be useful for future work on Majo- 
rana spin liquids, wc elaborate here on the differences be- 
tween our approach and that of BFLS. In particular, we 
discuss how fluctuations beyond mean-field theory may 
be incorporated in the approach of BFLS. We argue that 
following the standard route to construct an effective lat- 
tice gauge theory - a route that is successful in our for- 
malism - is problematic, except when symmetry allows 
for a partition of the lattice into dimers. Similar issues 
arise in constructing a projected wavefunction. However, 
we discuss a different route by which fluctuations may be 
included,^^'^'' which appears to be sound. This approach 
has, to our knowledge, not been studied in detail, and 
it will be very interesting to do so in future work. The 




FIG. 8. Orientation of nearest-neighbor bonds on the trian- 
gular lattice, used to define the mean-field Hamiltonians of 
Eqs. (D2) and (D5). This pattern of orientations respects the 
full translation symmetry of the triangular lattice. 

question of constructing a projected wavefunction, start- 
ing from the mean-field theory of the BFLS state, is still 
open. 

These differences notwithstanding, our approach can 
describe a very similar state to that of BFLS, differing 
only in the presence of a gapless s-fermion at the mean- 
field level. This state and the BFLS state are expected to 
have very similar physical properties. Indeed, it is con- 
ceivable that these two mean-field theories are different 
limits of the same phase. 

The starting point of Ref. 13 is a parton representation 
developed in Rcfs. 43 45, where the spin operators on 
each lattice site are represented directly using a triplet of 
S = 1 Majorana fermions: 

Si. = -yi%tl, (Dl) 

where {t\.,t'^,} = 26'^^5rr'- Focusing here and below on 

the triangular lattice and the BFLS state, these fermions 
are taken to obey a mean-field Hamiltonian 

{rr') 

where the orientations of nearest-neighbor bonds (rr') 
are chosen as shown in Fig. 8. As noted in Ref. 13, 
this mean-field Hamiltonian preserves the full trans- 
lation symmetry of the triangular lattice, but breaks 
time-reversal and certain point group symmetries. The 
fermions have an interesting locus of gapless excitations 
in momentum space. ^"^ In particular, the fermions are 
gapless at the F point of the Brillouin zone, and the 
vanishing gap at this point is protected by translation 
symmetry. 

As usual, Eq. Dl itself does not completely define the 
parton representation of the spin model; the physical 
Hilbert space must also be specified as a constraint on 

the larger parton Hilbert space. Because there arc an 
odd number of Majorana fermions on each lattice site, 
it is not possible to define a parton Hilbert space for 
individual lattice sites. Ref. 13 describes how one may 
proceed: the lattice can be (arbitrarily) partitioned into 
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dimers. Suppose that ri and r2 arc the two sites in some 
particular dimer. We can define three complex fermion 
operators by writing 

= (D3) 

which satisfy anticommutation relations 
{Crij-2' ("^rira)^} = ' ^e Can choosc the physical 
Hilbert space of the dimer to be those states with an 
even number of c*-fermions; the state with no fermions is 
the spin singlet, and the three states with two fermions 
form the spin triplet. This choice is enforced by the 
local Z2 constraint 

It should be noted that the choice -D^iT-a = — 1 is also 
legitimate; this picks out the states with an odd num- 
ber of c*-fermions, where the single-fermion states make 
up a spin triplet, and the unique three- fermion state 
is a spin singlet. These two choices of constraint are 
in fact interchanged by the Z2 "gauge transformation" 
tj.^ — >■ t^ii ~^ ~^T-2- This transformation does not 
satisfy the usual requirement for a gauge transformation 
that it leave the physical Hilbert space invariant. There is 
also a more conventional Z2 gauge freedom under trans- 
formations tj. -> TTrtl, where tt^ takes values ±1 and is 
constant on each dimer. 

For a given partition of the lattice into dimers, we can 
specify a sector of Hilbert space by choosing D^irg = ±1 
on each dimer. Importantly, because the spin operators 
are invariant under the Z2 transformations that change 
the constraint, the physics is the same in every sector. 
That is, all correlation functions of spin operators are 
the same in every sector. 

So far, we have described an exact parton representa- 
tion of the spin model, with no approximations. Now, we 
consider the mean-field starting point of Hq, and discuss 
how fluctuations may be included. We shall first follow 
the standard route to construct an effective lattice gauge 
theory, and show that this approach is problematic, ex- 
cept when symmetry allows for a natural partition of the 
lattice into dimers. 

We begin with a particular sector of the parton Hilbert 
space, with a fixed set of local constraints. For simplic- 
ity of discussion we choose Drir2 = 1 on each dimer. 
Those Z2 transformations preserving the constraint (i.e., 
those transformations that arc constant on dimers) play 
the role of gauge transformations. With this in mind, wo 
should define Z2 gauge fields ct^ and cr^, not on the links 
of the original lattice, but on links connecting nearby 
dimers. The local constraint -Drir2 = 1 is promoted to 
the gauge constraint Dt-it-2 =11^^! where the product 
is over those Z2 electric fields touching the dimer 
Hopping terms within a dimer are unchanged 
from Ho, but those between dimers are multiplied by the 
Z2 vector potential a^. Because any dimer covering of 
the triangular lattice breaks translation symmetry, it is 
clear that this effective gauge theory has less symmetry 



than Hq, and moreover the breaking of lattice symme- 
try depends on the dimer covering chosen. Therefore the 
mean-field Hamiltonian Ho does not properly capture the 
physics described by the more complete effective lattice 
gauge theory. Of course, if Hq itself breaks translation 
symmetry in such a way that there is a natural dimer 
covering, then we can employ this dimer covering in the 
construction of the effective gauge theory, and no prob- 
lems arise. 

The issue that arose here is essentially that, while the 
parton representation employed in Ref. 13 gives an ex- 
act representation of the spin model when the constraint 
-Drii-2 = 1 is imposed exactly, some symmetries are ex- 
plicitly broken the moment the constraint is softened, 
as occurs in the route we described to construct an ef- 
fective lattice gauge theory. A very similar issue arises 
and has been discussed in U(l) slave-rotor theories of 
the Hubbard model on bipartite lattices (see Sec. IV. A 
of Ref. 38). 

Similar issues arise in the construction of a projected 
wavefunction from the mean- field starting point of Ho- 
Given a partition of the lattice into dimers labeled by D, 
we define Vd to project onto the physical Hilbert space 
of each dimer. Then a projected wavefunction is given 
by IV') = T-'\''Po), where \ipo) is the ground state of Ho, 
and P = Ylo^D- Because P does not commute with 
translations, we expect that lip) is not invariant under 
translations and thus has less symmetry than Hq. Once 
again, if translations are broken in Hq so that there is a 
preferred dimer covering, this dimer covering can be used 
to construct P, and no issue arises. 

Despite these difficulties, fiuctuations can be included 
by a different route that does not suffer from the above 
issues. "'^'"'^ We can begin with the Heisenberg spin model, 
represented in terms of partons using Eq. (Dl), but with- 
out any local constraints. In the resulting Grassmann 
functional integral, wc have only the time-derivative term 
for the fermions, and the quartic spin-spin interaction. 
To obtain this description we have summed over all the 
sectors of Hilbert space described above; this is legiti- 
mate because the physics is the same in every sector. The 
spin-spin interaction can then be decoupled by standard 
means and Ho can be obtained as a mean-field saddle 
point. To our knowledge, fluctuations about this mean- 
field saddle point, or, indeed, about any saddle point in 
this construction, have not yet been studied. It will be 
important to investigate this in future work. Moreover, 
it may be possible to use these ideas to construct a pro- 
jected wavefunction for the BFLS state that does not 
suffer from the issues described above. 

Finally, we turn to the description of a state very simi- 
lar to the BFLS state, using our formalism. We consider 
the mean-field Hamiltonian 

Ho = ix* X] *^ ■ + X] 

{rr') (rr') 

where again the orientations of nearest-neighbor bonds 
(rr') are chosen as in Fig. 8. If it were possible to gap 
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out the s-fermions, without breaking more symmetries 
than are aheady broken in Hq, wc would obtain another 
description of the BFLS state. However, translation sym- 
metry requires both t- and s-fermions to be gapless at the 
r point, and this state is thus distinct from the BFLS 
state, at least at the mean-field level. It is conceivable 



that this mean-field state and the BFLS state could be 
two different frec-fcrmion limits of the same phase. Even 
if the two states are distinct beyond mean-field theory, 
this state is quite similar to the BFLS state, and we ex- 
pect that its physical properties are very similar to those 
elucidated in Ref. 13. 
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